Effects of the trapping potential on a superfluid atomic Fermi Gas 
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We examine a dilute two-component atomic Fermi gas 
trapped in a harmonic potential in the superfluid phase. For 
experimentally realistic parameters, the trapping potential is 
shown to have crucial influence on various properties of the 
gas. Using an effective hamiltonian, analytical results for the 
critical temperature, the temperature dependence of the su- 
perfluid gap, and the energy of the lowest collective modes are 
derived. These results are shown to agree well with numerical 
calculations. We furthermore discuss in more detail a previ- 
ous proposed method to experimentally observe the superfluid 
transition by looking at the collective mode spectrum. Our 
results are aimed at the present experimental effort to observe 
a superfluid phase transition in a trapped atomic Fermi gas. 
Pacs Numbers: 03.75.Fi, 05.30.Fk, 

With the increasingly impressive experimental ability 
to trap and cool atoms with fermionic statistics, a new 
and very active area within the field of ultracold atomic 
gases has emerged. Several groups are now involved in 
trapping and cooling atoms with Fermi statistics and 
temperatures as low as ~ 0.2Tp have been obtained for 
40 K and 6 Li [1]. One major motivation for this exper- 
imental effort is to observe a phase transition to a su- 
perfluid state predicted to occur below a certain critical 
temperature T c for a two-component gas [2] . 

In this paper, we examine various properties of such 
a two-component Fermi gas with an equal number of 
atoms N a of mass m in two hyperfine states \a =T>!) 
trapped in a spherically symmetric harmonic potential 
of the form Uo(r) — ma>f.r 2 /2. Assuming that the gas 
is in the dilute regime, the interaction between atoms in 
the two different hyperfine states can be well described 
by a contact potential gS(r) with g = Ai:h 2 a/m and a 
the s-wave scattering length. The dilute regime is de- 
fined as k F \a\ <C 1 with kp — ^/2mjxp/h and fip the 
chemical potential. For a < the effective atom-atom 
interaction is attractive and the gas is unstable toward 
the formation of Cooper pairs below a critical tempera- 
ture T c [2]. However, the specific formula for the criti- 
cal temperature was derived for a homogeneous system 
and is not necessarily valid for a trapped gas as we shall 
see. Indeed, it was shown in ref. [3,4] that the nature 
of the superfluid state depends on the coupling strength 
and the number of particles trapped. For many parti- 
cles trapped and/or stronger coupling strengths (keeping 
k F \a\ <C 1), the size of the Cooper pairs given by the co- 
herence length £,bcs = h 2 k F /(nmA), with A being the 
BCS pairing gap (for T = 0), is much smaller than the 
size R = y2/7p /m/ujT of the system. In this limit the 
system can be described as a quasi-homogeneous bulk 
superconductor with a locally varying chemical potential 



Mf(?") = fJ-F — Uo(r) [2]. For weak atom- atom interac- 
tions or for few particles trapped on the other hand, the 
size £,bcs of the Cooper pairs is comparable to the size 
of the system. The pairs are then formed between par- 
ticles with time-reversed angular momentum quantum 
numbers (l,m) and (l,—m) residing in the same har- 
monic oscillator shell and the gas is in the so-called in- 
trashell regime [3]. It can in this limit in many ways be 
regarded as a large nucleus with some properties quali- 
tatively different from the bulk regime. The transition 
between the two different regimes is roughly determined 
by ^bcs — R or equivalently A ~ Hlot- Using the bulk 
theory prediction A = (2/e) 7 / 3 ^ exp(— Tr/2k F \a\) [5] 
and the Thomas-Fermi result \x F = (6A r (J ) 1 / 3 fiwr, the 
equation A = Tilot yields N* — g (|) 7 exp(37r/2fci?|a|). 
For N a <C N* atoms trapped, the gas is in the in- 
trashell regime and for N a ^> N* atoms trapped it 
is in the bulk regime. With k F \a\ = 0.1 this condi- 
tion yields N; = 4 x 10 20 and for k F \a\ = 0.3 we get 
N* = 9 x 10 6 . Current experiments have N ~ 0(1O 5 ) 
particles trapped [1] suggesting that in the dilute limit 
one is likely to be in the intrashell regime. The ex- 
periments seem to be aimed at using optical traps as 
the confining potential since these offer more flexibility 
with respect to which atoms can be trapped [1]. These 
optical traps are rather tight with trapping frequencies 
lu t ~ O(10 3 Hz) yielding huj T ~ O{10 2 nK). Thus, the 
condition k B T c < frw T [k B T c < k B T F = (GTV^) 1 / 3 ^] 
for the intrashell regime is less restrictive than for mag- 
netic traps. However, cooling schemes as the one de- 
scribed in [6], which enables the lowering of T/Tp by 
orders of magnitude, will probably have to be employed. 

For kp\a\ > 0.3 the gas is no longer in the di- 
lute regime. Recently, a number of papers have been 
published examining the gas in this strongly correlated 
limit [7], It is presently unclear however, whether the gas 
is stable toward spinodal decomposition or other insta- 
bilities and it is outside the scope of the present paper to 
describe the gas in this regime. 

One property of great experimental importance is the 
critical temperature T c which we will now estimate for the 
gas in the intrashell regime. In this limit there is only 
pairing between particles with quantum numbers (n, I, m) 
and (n, /, — m) with the unperturbed harmonic oscillator 
shell energy e„ = (n+3/2)fkuT- This fact, combined with 
a number of semiclassical approximations allows one to 
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derive an effective Hamiltonian of the form [3] 



2FlW T G i\m+m'»t at 



nlrn ,n' V m' 



(i) 



where f2„ F = (np + l)(np + 2)/2 is the unperturbed 
degeneracy of the highest harmonic oscillator shell occu- 
pied with energy e„ F (ignoring the Hartree field for the 
moment), a n i m is the operator removing a particle in the 
shell with energy e„ and angular momentum quantum 
numbers l,m and £„ = e n — fj,p. The effective coupling 
strength giving rise to the pairing correlations is given by 
G = 327T - 2 A;^ |a| /15. From Eq. (1) a number of analytical 
results can be derived. The gap equation becomes [3]: 



G 



2f(E n )] 



(2) 



where we have as a first approximation to the pseudopo- 
tential scheme introduced the cutoff n < 2np. The 
quasiparticle energy is E n = (£ 2 + A 2 ) 1 / 2 and f(x) = 
[exp(f3x) + with (3 = l/k B T. As A « frw T in the 
intrashell regime, shell effects are pronounced; the be- 
haviour of the system depends crucially on whether the 
system is in a filled shell configuration or whether the 
highest occupied shell is only partially filled [3] . Since the 
filled shell configuration occurs only for certain "magic" 
numbers of particles trapped, we expect the typical ex- 
perimental situation to be that of a partially filled high- 
est shell [np — (np + 3/2)Hujt]- In this case, and for 
ksTc <C tux>T appropriate for the intrashell regime, Eq. 
(2) yields the critical temperature 



G 



Ao 



huj T 2 - Un(e^n F )G 2frtu T ' 



(3) 



Here 7 ~ 0.577 is Euler's constant and Ao is the gap for 
T = in the intrashell regime [3]. For G — > 0, we have 
ksTc/hujT — G/2. In this very weak interaction limit, 
there is only pairing in the partially filled shell right at 
the chemical potential (n — np) and the shells away from 
the chemical potential are unaffected by the pairing. For 
increasing G, the term ln(e 7 nF)G in Eq. (3) becomes 
important and kBT c /hujT > G/2. Physically, this corre- 
sponds to the interaction now being so strong that there 
is significant pairing in the shells away from the chemical 
potential; each Cooper pair is however still only formed 
between particles within the same shell. We see that T c 
in the intrashell regime given by Eq. (3) is very different 
from ksT c oc fip exp(— ir/2kp \a\) relevant for the bulk 
regime [2,5]; for G -C 1, Eq. (3) yields a critical temper- 
ature which is orders of magnitude higher than the bulk 
result. This stems from the fact that the shell structure 
of the trap levels increases the density of states for certain 
energies. For ln[exp(7)n^]G = 1/2, Eq. (3) diverges but 



this is clearly unphysical as the intrashell limit is charac- 
terized by ksT c <C Hlot- For such large values of G the 
intrashell pairing ansatz is invalid and Cooper pairs are 
formed between particles residing in different shells. It 
should be noted that Eq. (3) does not include the effect 
of induced interactions coming from density oscillations 
of the gas. In the homogeneous case, this effect is known 
to reduce T c by a factor ~ 2.2 [5] and it should be in- 
vestigated which influence this effect has in the intrashell 
regime as well. However, it must be expected that the 
critical temperature including the induced interactions 
will be of the same order of magnitude as Eq. (3) and 
thus very different from the bulk theory prediction. 

To check the validity of Eq. (3), we plot in Fig. 1 
kBT c /Hu>T as a function of G. The solid curve is Eq. (3) 
and the x 's are obtained from a numerical calculation 
solving the Bogoliubov-de Gennes (BdG) equations de- 
rived from the full Hamiltonian using a pseudopotential 
scheme [8]. The dashed line depicts the bulk prediction 
ksT c oc fip exp(— ir/2kp\a\). We have chosen np ~ 30 
giving N = 2N a ~ 10 4 particles trapped. There is excel- 
lent agreement between the numerical results and Eq. (3) 
for G < 0.05. We conclude in agreement with ref. [3] that 
the intrashell description works well for ksT c <C hwp. In 
this limit on the other hand, the bulk theory is qualita- 
tively wrong as expected. For stronger coupling, the na- 
ture of the superfluid state undergoes a transition and the 
intrashell description breaks down. There is now pairing 
between particles residing in different harmonic oscillator 
shells. This breakdown of the intrashell theory occurring 
for ksT c > 0(Hujt) for larger coupling strengths G is 
not depicted in Fig. 1. Eventually the coherence length 
S,bcs becomes much smaller than the size of the system. 
As explained above, for N ~ 10 4 trapped the gas is no 
longer in the dilute limit when £,bcs ^ R- 
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FIG. 1. The critical temperature kBT c /hu)T as a function 
of G. +'s and x's are with and without the Hartree field 
respectively. 

The Hartree field introduces a dispersion on the Har- 
monic oscillator levels with respect to the angular mo- 
mentum / lifting the shell degeneracy fl n = (n + l)(n + 
2)/2. When np < 30 (depending on the size of the cou- 
pling strength), this splitting is small compared to the 
pairing energy A and it can be ignored [3]. For such 
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systems, Eq. (3) for k B T c should be accurate. For very 
large systems with rip ^> 200 on the other hand, the 
Hartree splitting of the shell is much larger than A and 
there is for weak coupling only pairing within states with 
a fixed n and I [3]. In order to examine the effect of 
the Hartree field, we plot in Fig. 1 as +'s k B T c /hu>T as 
a function of G obtained from a numerical calculation 
where the Hartree field is included. We see that effect of 
the Hartree field is to suppress T c for small G due to the 
fact that the lifting of the shell degeneracy reduces the 
density of states at the chemical potential. For stronger 
coupling, the Hartree field increases T c as it in this limit 
enhances the spatial overlap between the orbitals of the 
Cooper pairs. However, from Fig. 1 we see that Eq. (3) 
agrees qualitatively with the numerical results of the crit- 
ical temperature. This is because uf ~ 30 is not quite 
in the limit where the Hartree field can be neglected, but 
also far away from the regime uf ^ 200 where there is 
only pairing within levels with a single lvalue. Again, 
these results are in agreement with earlier calculations 
on the T = properties of the gas [3] . 

We now examine the temperature dependence of the 
pairing energy A in the intrashell regime. Solving Eq. 
(2) for A(T)/huT < 1, we obtain 



T 



A(T) 



k B T c In 



l+A(T)/2fc B r c 



l-A(T)/2k B T a 



(4) 



from which A(T) can be found by inversion. For (T c 
T)/T c < 1 such that A(T)/k B T c < 1 we have 



A(T) 

k B T c 



= VlVl - T/T c . 



(5) 



In Fig. 2, we plot the lowest quasiparticle energy 
E{T)/K B T C = A(T)/k B T c (since £„ F ; = 0) as a func- 
tion of T/T c . The solid line is Eq. (4) and the x's and 
+'s are numerical results from a solution of the BdG Eq. 
with ~ 10 4 particles trapped and G = 0.067 where the 
effect of the Hartree field has been excluded for the x 's. 
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FIG. 2. The gap A(T)/k B T c as a function of T/T c . +'s 
and x 's are with and without the Hartree field respectively. 

We see that there is good agreement between the nu- 
merical results and Eq. (4) in the intrashell regime. In 
particular, we have A(T = 0)/k B T c ~ 2 as was predicted 



in Eq. (3). Again, the effect of the Hartree field does not 
change the qualitative behaviour of the system for this 
set of parameters since rip <C 200. So we conclude that 
Eq. (4) describes A(T) well in the intrashell regime. 

An important question is how to experimentally ob- 
serve the superfluid transition. One way is to look at the 
collective mode spectrum of the gas on which superflu- 
idity should have several observable effects [9,10]. When 
the gas is in the intrashell regime, two such effects are 
the emergence of low energy density oscillation modes 
of monopole and quadrupole symmetry. The energy of 
these new modes should scale as Huj ~ 2 A <C Hujt [9] . To 
examine this, we plot in Fig. 3 as +'s the frequency ojq of 
the low energy quadrupole mode as a function of T. The 
collective response is calculated using a RPA-scheme in- 
cluding the superfluid correlations as described in ref. [9] . 
The solid line gives 2A(T) obtained from Eq. (4) and we 
also for comparison plot as x's 2A(T) obtained from a 
numerical solution to the BdG eqn. The parameters used 
are the same as for Fig. 2 and the effect of the Hartree 
field is included. 
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FIG. 3. The low energy quadrupole frequency huiQ/kBT c 
as a function of T/T c . The +'s are obtained from a numerical 
RPA calculation, the x's are 2A/hu>T from a numerical solu- 
tion of the BdG eqn., and the line is 2A/!ujt with A given 
by Eq. (4). The inset shows S(u) (in trap units) as a function 
of oj/ut for T = 0. 

We see that there is fairly good agreement between the 
numerical results from the RPA calculation and Eq. (4). 
The energy of the collective mode given within the RPA 
is approximately 2A indicating that the mode consists of 
breaking a Cooper pair. The fact that the collective re- 
sponse frequency in general is lower than the pair break- 
ing prediction 2A can be understood as follows: The 
excitation consists of the breaking of a Cooper pair of 
monopole symmetry costing 2A in energy. However, the 
broken pair of particles then form a new Cooper pair 
with quadrupole symmetry; the difference 2A — hujQ can 
be interpreted as the binding energy of the quadrupole 
Cooper pair. The excitation can therefore be understood 
as the breaking of a Cooper pair with monopole sym- 
metry and the subsequent formation of a Cooper pair of 
quadrupole symmetry, i.e. a pair vibration mode. 

The frequency ljm of the lowest monopole mode is ap- 
proximately the same as for the lowest quadrupole mode 
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shown in Fig. 3. Again, the excitation frequency is lower 
than the simple pair breaking prediction 2A because 
the broken pair form a new Cooper pair in an excited 
monopole state. We will in a future publication examine 
in more detail the binding energy of these excited Cooper 
pairs of monopole and quadrupolc symmetry giving rise 
to the collective mode frequencies being somewhat lower 
than 2A. For the purpose of the present paper however, 
it is adequate to conclude that the onset of superfluidity 
below T c indeed does give rise to new density oscillation 
modes of monopole and quadrupole symmetry with en- 
ergy Hlo ~ 2 A in the intrashell regime. 

The collective modes can be excited by modulating 
the trapping frequency producing a perturbation with 
the right symmetry [11]. There is however one im- 
portant point concerning the monopole pair vibration 
mode: In order to excite this mode, one should ap- 
ply a perturbation of the trapping potential 5H = 
sin(wi)E CT / d 3 rF(r)V>t( r )^ (T ( r ) w ith F(r) ^ ar 2 . This 
is because in the intrashell limit, the Cooper pairs are 
only formed between states within the same harmonic 
oscillator shell with unperturbed energy (n + 3/2)fkdT- 
Neglecting the effect of the Hartree field, the states within 
one shell all have the same expectation value for r 2 . 
Thus, due to the orthogonality between the eigenstates 
of H, a perturbation oc r 2 (such as the modulation of 
the trapping frequency) will not probe states which cor- 
respond to excitations of pairs within the same harmonic 
oscillator shell and the monopole pair vibration mode 
fkJM ~ 2A will not be excited. To excite this mode, the 
scheme using off-resonant lasers demonstrated might be 
useful [12]. For the numerical calculations, the monopole 
mode was excited by a perturbation with F(r) oc r 4 . 
A related issue is the spectral weight with which the 
pair vibration modes are excited as compared to the 
normal phase modes at u <~ 2ujp- The spectral weight 
#M « EJ(0|&ff|0)| 2 $(fiw - -Bo - E v ) with |0) and \v) 
being the ground and excited states with energies Eq and 
E v respectively, is for T = in the quasiparticle approx- 
imation 

S(uj) oc | / d 3 r[u v v v i + v ri u r , / ]F\' 2 d(hu) - E v - E v >), 
vv' (6) 

where the elementary quasiparticles with excitation en- 
ergies E v are described by the Bogoliubov wave func- 
tions u v (r) and v v (r) [13]. r\ = (n,l,m) is shorthand for 
the quantum numbers characterizing the quasiparticle. 
The collective pair vibration is predominantly built out 
of quasiparticle excitations with n = n' = np whereas 
the normal phase excitation u <~ 2lot is built out of ex- 
citations with n = np — 2, np — 1, np and n' = np + 2 
or vice versa in Eq. (6). Now, for the n = np shell, 
even a minute pairing will mix the particles and holes 
yielding u nF i m (r) = v nF i m (r) = 2- 1 / 2 ^ F , m (r) where 
<f>n F im( r ) are t ne non- interacting harmonic trap eigen- 
functions with energy (n + 3/2)Hwt [3]. This makes 
the UrjVrj' in Eq. (6) large and depending on the spa- 



tial overlap between F(r) and the excitations involved, 
the strength of the pairing vibration will be of the same 
order as the normal phase excitation to ~ 2u>t, although 
in general somewhat smaller since it is built from excita- 
tions from only the n = np shell whereas the lo <~ 2iop 
mode is built from the shells rif — 2,np — l,np. This 
is illustrated in the inset in Fig. 3, where we plot S(u) 
(calculated using the RPA) for T = with the param- 
eters chosen above for the quadrupole mode excited by 
F(r) = r 2 Y2o{0, <p). We see that the strength of the pair- 
ing vibration is comparable to the normal phase mode 
thereby making it experimentally observable. 

When the trap is not spherical, the single particle en- 
ergies and wave functions are different making the quan- 
titative predictions of this paper invalid. However, if 
all trapping frequencies iOi (i — x,y,z) are such that 
hwi 3> ksT c , the qualitative conclusions of this paper 
still holds: T c can be orders of magnitude higher than 
the bulk result due to increased density of states for 
certain energies and there will be pair vibrations with 
fku ~ 2A. If, on the other hand, one frequency is loose 
(say tko z <C ksTc), then the gas will be in some interest- 
ing mixture of the bulk and the intrashell regime which 
warrants further study. 

In conclusion, we have shown that for experimentally 
realistic parameters, the trapping potential has a signif- 
icant influence on the properties of a trapped atomic 
Fermi gas in the superfluid phase in the dilute limit. 
In the intrashell regime, we used an effective Hamilto- 
nian to derive analytical results concerning the critical 
temperature, the temperature dependence of the pair- 
ing gap and the lowest collective oscillations of monopole 
and quadrupole symmetry. These results were shown 
to agree well with numerical calculations. We also dis- 
cussed in further detail a previously proposed scheme to 
detect the onset of superfluidity with decreasing temper- 
ature by looking at the collective mode spectrum. Much 
of the analysis of this paper should be relevant for the 
present experimental effort aimed at observing the su- 
perfluid transition in a two-component trapped atomic 
Fermi gas. 

We acknowledge valuable discussions with C. Salomon, 
B. Mottelson, and J. R. Krumrinc. 
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